Characteristic features in non-Markovian noise spectrum of transport current 
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Based on the construction of an efficient scheme for non-Markovian noise calculation, we analyze 
the noise spectrum of transport current through interacting single quantum dot and double dots. 
We show that the non-Markovian correction is remarkable, ft leads to a number of characteristic 
spectral structures, including that due to finite bandwidth and that sensitive to and enhanced by 
the magnitude of Coulomb interaction. 
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Shot noise due to charge discreteness in mesoscopic 
transport has stimulated great interest in recent years. 
It provides much rich information beyond the average 
current p], Conventionally, shot noise and higher 
cumulants of current in full counting statistics (FCS) 
are largely restricted to zero frequency, and the Born- 
Markov memoryless master equation approach is em- 
ployed @, i, i, i, h, s s m tu m m q m m. 

Memory effects on the first few cumulants of FCS at 
zero frequency were investigated recently, revealing that 
the non-Markovian correction becomes increasingly im- 
portant to higher cumulants [l7j . Non-Markovian effect 
should be more pronounced at finite frequency, especially 
in the high-frequency or quantum regime. 

Memory effect manifest primarily in the short time dy- 
namics. It supports some energy-unconserved transient 
processes as the energy uncertainty is significant there. 
The relevant time scale is however physically related to 
the electronic structure of the many-particle memoscopic 
system and dynamical property of system-reservoir cou- 
pling statistics such as finite bandwidth. Some trans- 
port dark channels in steady state current and/or zero- 
frequency current cumulants related to long-time dy- 
namics will be activated. As results, the finite Coulomb 
interaction and finite bandwidth will manifest their sig- 
natures in the noise spectrum. 

In this work, instead of the infinite Coulomb interac- 
tion and wideband limit ( WBL) [l8|, H(| , we demon- 
strate some remarkable features resulted from the in- 
terplay between the non-Markovian memory effect and 
intrinsic electronic structure of system-leads complex. 
Upon the establishment of an efficient scheme for non- 
Markovian noise calculation, we investigate the circuit 
quantum noise in the full frequency range, covering over 
the Coulomb interaction and even finite bandwidth. A 
number of striking features will emerge. For example, 
provided that the energy-allowed transport channels have 
been defined, the Markovian noise spectrum no longer re- 
flects the magnitude of Coulomb interaction [2l[ , but the 
non-Markovian spectrum does contain that information 
and even is enhanced by it. The finite bandwidth will also 
lead to qualitatively different result from the Markovian 
treatment. These prominent features provide a very use- 
ful access in experiments to probe important parameters 



of devices. In addition, they are also important to im- 
prove the signal-to-noise ratio in the regime where the 
fluctuation energy is larger than the charging energy for 
the electron transport setup as readout devices for solid- 
state charge qubits [201 ] - 

Our study is based on the particle- number "n"- 
resolved quantum master equation (QME) approach 
8, 14], implemented with the memory kernel prescrip- 
tion. Let us start with the conventional QME for the 
reduced system operator containing memory [22| : 



p(t) = -iCp(t) 



drE(t, t)p(t), 



1) 



with £(t,r) = £(*-t) = <£ , (i)e- l£ ( t - r )/: , (T)) B . Here, 
C ■ = [H, ■ ] is the system Liouvillian; C'(t) ■ = [H'(t), ■ ] is 
the system-reservoirs tunnel coupling Liouvillian; (■ ■ -) B 
denotes an average over the reservoirs bath degree of free- 
dom of electrodes with different chemical potentials un- 
der bias voltage. The electron reservoirs bath Hamilto- 



nian assumes h B = e a ki 



a a k 



The system H is rather 



general and written in term of the creation {d^ = d^} 
and annihilation {g? m = d~} operators of electron in sys- 
tem obitals. The tunnel coupling in the /i B -interaction 
picture is H'(t) = E^/'^'tWaL^ + H^)^* 8 ' = 

E a) i [/ij' + dlfi^\t)], in the view point of dissi- 
pative dynamics. The effect of leads on central system 
is completely characterized by the reservoir correlation 
functions C$ v {t - r) = (/^ (t)) b , where a takes 
the opposite sign of a = + or — . It results in the self- 
energy kernel E(f) in Eq. JTJ) of 



nt)6 = Y,{l^ D{ 3& c )°\ + H - c -} 



(2) 



where £) = £„ e- iCt C^(t)dZ. The above QME 

invokes the second-order Born approximation and is 
valid for sequential tunneling. For later use, let us in- 
troduce the half-Fourier transform that resolves Eq. (JTJ) 
with the formal solution p(u) = [Ti(u>)+i(jZ — ui)] p(to). 
The stationary state is then p = [—iu}p(ui)] u =o- The res- 
olution of Eq. f2| reads 

t(u)6 = ]T [d%D$(£ - oj)6 - OD^(C + <jj\ , (3) 

a/j.er 



2 



with D$[(u) = £„ J^dte-^C^MtK- 

Rather than the above unconditional QME, a richer 
information contained equation for conditional state will 
be more desirable. This is the well established particle- 
number-resolved QME, which contains such transport in- 
formation as current, shot noise, and even all the higher 
moments of current fluctuations [Til ]. The key quantity 
in the "n" -resolved QME is p( n "\ the reduced system 
state conditioned on n a electrons passed through the 
tunnel junction between the a-th lead and the central 
system. The unconditional state is related to p( n °) via 
P = E„ Q P ina) - Following Refs. @, HI [H, especially 
the idea of decomposition of the Hilbert space presented 
in 14], the non-Markovian version can be readily formu- 



lated out as (setting = J2 a ^ol* an< ^ a ' ^ 



(t-r;C) 



xp^\r)~Dl%(t-r-C)p^\r)dl 



D^(t-r;£)p 



(n a +al) 



(r)d° 



H.c 



}• (4) 



The conditional state is straightforwardly related to the 
distribution function for the tunnelled electron num- 
bers: P{n a ,t) = Tr[p("°)(i)]. With the knowledge of 
P(n a ,t), all transport properties can be obtained. First, 
for the current, through I a (t) = —■^^ na n a P{n a ,t), 

we have I a (t) - - ft ^t(-)i 



ja \c,uj)p(uj), with 



■ f ta dTjy>{C,t- t) P {t) or J a (w) 



d^dD^iC + u)] ±dl [dD(+)t(£ + c)]}. (5) 



+ 



The stationary current reads I a = —J& \c,u> = 0)p. 

For the circuit current noise spectrum, we have S(lu) — 
oS'l(w) + bSji(u>) — abS c (to). The capacitive coupling 
related parameters a and b link the total circuit cur- 
rent with the left and right junction currents as I(t) = 
alhit) — bln(t). Accordingly, Sl/r(w) is the noise 
spectrum of lL/n{t), and S c (ui) is the charge fluctua- 
tion spectrum on the central dots, satisfying Q(t) — 
~[Ih(t) + Init)]- For Sl/rM, it follows the MacDon- 
ald's formula S a (u>) — 2u; dtsmcjt-^(n'^(t)) , where 
(ni(t))=Z na nlTrp( n ^(t)=i: na niP(n a ,t). With the 
help of Eq. ((4J and the MacDonald's formula, the cur- 
rent noise spectrum can be finally expressed as 



(6) 



5 Q (w) =4wIm{Tr[ l 7 t (-)(£,a;)^ a (a;)]} 
+ 2Re{Tr[jM(£,u)p}}, 

where N a (uj) = i[£(u>) + i{£ - u))]^ jt ] {£, u)p/u. 

For the charge noise, the current conservation leads 
to S c (lo) = 2S L rM + S l (uj) + 5r(£j), with Slr(w) = 



2ulm{Ti[J^>{£,u;)N R (to) + J^~'{£, lu)N l {lu)}} re- 
sulted from the MacDonald's formula as described above. 

To illustrate the non-Markovian reservoir coupling ef- 
fect on noise spectrum, we choose a simplified Lorentzian 
spectral density, T af , u {u;) = T a W%/[(uj - p a ) 2 + W%\, 
if it is nonzero. The WBL is achieved when W a — > 
oo. The reservoir correlation functions determine the 
rates of electron tunneling between the system and 
leads. They are related to the spectral density via 
the fluctuation-dissipation theorem (24|: Cailv(t) = 
/ duie aiut T atlv (bj)/[l + e ^(«-A»a)] ) w i t h p a the chemical 
potential of the a-lead and (3 the inverse temperature. 

Before going to applications, two general remarks are 
worthwhile, (i) The key feature of non-Markovian mas- 
ter equation is the memory kernel in Eq. |1| , which leads 
to the frequency-dependent tunneling rates. As an il- 
lustration, take a look at the term D^J{e^ — u))pd^ in 
Eq. ([3]). It describes the tunneling from the a-lead to 
the system state of energy e M , associated with an en- 
ergy of w absorbed (setting H = 1). On the other hand, 
the Markovian approximation is of w-independent tun- 
neling rates [2l| . (ii) The frequency-dependent tunneling 
rates have profound effect on the noise spectrum. For 
|w — Pa\ 3> W, the tunneling rates become zero, lead- 
ing thus to vanishing noise at high frequency limit. In 
contrast, the Markovian approach cannot predict this be- 
havior [2l|, as the tunneling rate exhibits no frequency 
dependent. Even for WBL model where both the non- 
Markovian and Markovian treatments do not approach 
to the zero noise result, qualitative differences arise, as 
will see below. 

Preliminary Insight. Consider the simplest system, a 
single dot with single spinless level, H = eoa^a, involving 
only two states, the empty |0) and the occupied |1). We 
also assume large bias and low temperature: /il > £o > 
Pb, and the Fermi function / Q (w) = 1 for u < p a , and 
zero otherwise. For the charge occupation fluctuation on 
the dot, S n (lu) = dt e iut (N(t)N(0) + N(0)N(t)), we 
find that (denoting 7 = Tl/Fr) 




< U) < UJlo 
7/2, UJ LQ < LO < CJRO 

7)/2, uj > cjro 



(7) 



where uj a0 = \p a - e |, and Sj}(u) = 4I/(T 2 + J 1 ) is 
the result under Markovian approximation, with T = 
Tl + Tr. Therefore, in frequency domain, the non- 
Markovian effect is manifested as a discontinuity at fre- 
quencies a>Lo(Ro)iJL e fl ec t m g the nature of resonance with 
the Fermi level [18[ . Figure Q] depicts the time domain 
counterpart of Eq. ([7]), i.e., the charge number correlation 
function Sjv(i) where nonexponential decay behavior is 
evident. Note that Sjv(i) is nothing but the probability 
for an electron remains in the dot at time t, assuming it 
was at t — 0. The nonexponential decay, which is typical 
for non-Markovian dynamics, essentially arises from the 
nonlocal time memory effect in Eq. Throughout this 
work, we set pi, = eV/2 and ^r = —eV/2. 

Interacting Single Dot. The relevant Hamiltonian 
reads H = e^n^+Un^hi, where — a^a^ and p the 
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FIG. 1: Charge correlation function Sjv(i), under non- 
Markovian (solid line) and Markovian (dashed line) treat- 
ments. Parameters (in unit of meV) are T = 0.05, eV = 5, 
e = 1, T = 1, and W = 15. 



spin index, so that C^},(oj) = C^(cj)5 ijv . To display 
the effect of Coulomb interaction more transparently, we 
assume the dot level spin-degenerate, ej = e^ = eo- Also, 
we consider the transport in strong Coulomb blockade 
regime where eo + U > /Xl > eo > Ma- 

Let us start with region of lo < u> a o and carry out the 
noise spectra analytically, which in terms of Fano factor 
S(w)/(2I) read (denoting r cff = 2r L + T R ) 



2(4r 



r 2 

1 R 



2r off n L 



2(4r 



r c 2 ff 

2 
L 



Ir4 



r R (rf ff + w >' 
) _ r cff 
^ 2 ' r L (r 2 ff + w 2 )^ 
x [n R + (lu 2 + ^ 2 r L /r cff - 6r 2 )$ R ] , 
2^ 2 r cff $[r ofr $ - r R w] 



r 2 

1 off 



(8a) 



(8b) 



r 2 

1 off 



r L r R (r 2 ff + ^ 2 ) 



where U a = w$ Q $ + 2r L r R # + 2r 2 ,* Q + uj 2 T a ,<$> a /T cS , 
$ Q (e ,cj) = x«(eo-^)-Xa(eo+^) = $a, * 



Ma)] (*iS 



$ = $ L -$ R , and Xa ( a; ) = ^Re*[I + ^(x 
the digamma function). At zero frequency, $ a (eo, 0) = 
the second terms in Eq. ([5]) vanish, and the Markovian 
result (the first term) is recovered [21[ . The second terms 
in Eq. ([8]) , which arise in fact from the renormalization of 
the dot level, reduce the tunneling rates and suppress the 
shot noise. They can be neglected in the high-frequency 
regime considered hereafter. 

Increasing frequency usually results in non-Markovian 
effect more distinct, as the short timescale dynamics be- 
comes more important. Consider the noise spectrum at 
intermediate frequencies. From Fig. 2, we notice a num- 
ber of characteristic frequencies, i.e., ix>l(R)o = IMl(R) - e o| 
and wl(r)i = eo + U — /l«l( R ), at which sudden change of 
fluctuation magnitudes takes place. The processes in- 
volve electron transfer between dot and leads, associated 
with the energy absorption/emission of detection. The 
resonant step-structures are completely absent in Marko- 
vian treatment; much richer non-Markovian transient dy- 
namics occur in time domain. 

Markovian theory concludes also that in Coulomb 
blockade regime the noise spectrum does not depend on 




FIG. 2: Noise spectrum of transport current through an in- 
teracting quantum dot, decomposed to (a) junction currents 
and (b) charge-number fluctuation components, respectively. 
Parameters (in unit of meV) are U = 10, T = 0.05, eV = 6, 
eo = 1, r = 0.4, with finite bandwidth of W = 15. The insets 
are the corresponding results of WBL. 



U [2l|. This result seems plausible, since noise is de- 
termined by the random occupation of the dot, which 
would be no longer relevant to the magnitude of U, pro- 
vided the dot has already been in the strong Coulomb 
blockade regime. This argument is, however, valid only 
for long-time dynamics. With the increase of frequency, 
the energy nonconserved states will be involved, due to 
the energy uncertainty in the relevant short-time regime, 
and the noise becomes [/-dependent. 

In the WBL where the memory arises only from Fermi 
distribution, we obtain (i) F a (uo a o < ui < ui a i) = |(1 + 
7 Q ) and F a (u> > oj al ) = \ + j a , where 7l = 1/tr = 7! 



and (ii) F c = 1 + |( 7 + 7- 1 ), 1 + 7 + ±7^ and 1(3 + 

2 7 + f° r W R0 < u < ^li, luli < lu < w R i, and 

ui > w R i, respectively. These results are plotted in the 
insets of Fig.O The noise can be either super-Poissonian 
or sub-Poissonian, depending on the ratio 7 = Fl/F r . 
Strikingly, it differs from the Makrovian results: F^f = 
F^ = 1 (i.e., Poissonian noise) and F^ 1 = 2. 

Remarkable non-Markovian effect is also manifested at 
high frequency limit. For the finite bandwidth model, 
the non-Markovian theory predicts a correct vanishing 
noise, while the Markovian treatment cannot give this 
prediction. For infinite wide-band model, the contradic- 
tion is not so sharp, however, interesting difference also 
emerges. Consider, for instance, Fl and F R . The Marko- 
vian treatment gives i*L = -Fr = 1, i.e., the Poissonian 
behavior. However, the non-Markovian treatment results 
in Fl = 5+7 and F R = 1 + \^~ x for strong Coulomb in- 
teraction eo + U > fih > eo > /U R ; and i*L(R) = ^(l + T* 1 ) 
for weak Coulomb interaction /iL > eo,eo+U > /i R . They 
can be either super-Poissonian or sub-Poissonian depend- 
ing on the ratio 7. Finally, the results also show that the 
noise is enhanced by Coulomb interaction strength. 

Interacting Double Dots. Now consider the trans- 
port through two coupled quantum dots, described by 

H = E Q =L,R eQ °£ d ° 



-UN L N R +n(dUn + dtd L ). We as- 
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FIG. 3: The circuit noise spectrum of coupled quantum 
dots under non-Markovian (solid line) and Markovian (dashed 
lines) treatments. Parameters (in unit of meV) are the same 
as Fig.0 together with el = £r = £o = 1 and Q — 2.5. The 
inset is the corresponding result of WBL. 

sume that the intra-dot Coulomb interaction is infinitely 
large; thus only the inter-dot interaction term appears 
explicitly. Moreover, for simplicity, we consider spinlcss 
electron here. In the strong inter-dot Coulomb block- 
age regime, the charge configurations relevant to steady- 
state transport are the empty |0), left-dot occupied \L), 
and right-dot occupied \R) states. Similar to the single- 
dot case, the non-Markovian memory effect will also in- 
corporate the double-dot occupied \d) state into dynam- 
ics. The characteristic frequencies in noise spectrum are 
now determined by the respective positions of the eigen- 
encrgics E e and E g of the coupled dots from the Fermi 
levels, i.e., \E e , g — ji a \ and \E e ^ g + U — jj, a \. In addi- 
tion, another characteristic frequency, the Rabi frequency 
A eg = E e - E g = y/e 2 + 4fl 2 , with e = e L - (r, will 



strongly affect the noise spectrum, particularly for a large 
n. The resultant noise spectrum is depicted in Fig.[3J 

For the double dots, interesting non-Markovian fea- 
tures also appear at high frequency limit. For the WBL, 
the Fano factor reads 

F a {io) =r a [r| + 4(2 + 7-*) n 2 + 4e 2 rcn }/(8T R n 2 ), (9) 

where e rcn = £ + AE L - AE R , with AE a = \ [xa(eo) - 
X«( e o + U)], accounts for the renormalization effect. 
Strikingly, it follows Eq. that the noise can be highly 
super-Poissonian, if the coherent coupling f2 between the 
dots is small. This result can be understood as follows. 
At high frequency limit, S a — T a , while the steady-state 
current reads / = 4F R ft 2 /[r2 + 4(2 + j-^fl 2 + 4e 2 en ] . 
In other words, the noises of the left and right junction 
currents do not vanish, but the steady-state current will 
if the coupling O — ^ 0. Differing remarkably, the Marko- 
vian treatment results in — 21, thus = 1. 

In summary, upon the establishment of the "n"- 
resolved QME with memory, we illustrated that a non- 
Markovian treatment is essential to correctly account 
for the frequency-dependent shot noise. Markovian ap- 
proach is applicable only in the classical regime. It breaks 
down completely, even for the wide-band leads, as long 
as the measurement frequency becomes comparable to 
or higher than the characteristic excitation frequencies. 
This is the quantum regime where memory effect is strik- 
ingly prominent. The resultant characteristic features in 
the noise spectrum provide very useful access to some im- 
portant system parameters of study. Support from GRC 
Hong Kong and NNSF China is acknowledged. 
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